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Abstract 

The spin dependent deuteron structure functions are investigated in the rel- 
ativistic model of deuteron. It is shown that the deuteron structure function 
g^iX) can be considered properly only in the infinite momentum frame where 
the virtual photon has pure transverse components (q = —q 3 — ► at P — ► oo), 
whereas in the conventional Breit frame the impulse approximation for g^iX) may 
be violated due to qq pair creation by the virtual photon. It is shown that due 
to relativistic effects the deuteron structure function g^iX) has a nontrivial con- 
tribution determined by first derivative of the nucleon structure function gi(x). 
Such contribution is small but it increases when X — ► 1. It may achieve about 
10 % at X ~ 0.75 and could be essentially larger for the transverse structure 
function gj_(X). 

1 Introduction 

The investigation of the spin-dependent nucleon structure functions (SF) has attracted 
much attention recently. Experiments are prepared at SLAC(E142/143) [1], CERN(SMS) 
[2] and DESY(HERMES) [3]. The neutron SF can be extracted from the light nuclei 
( 2 D or 3 He) data, and nuclear effects should be taken into account to obtain accurate 
results. The procedure involves convolution of the nucleon infinite momentum frame 
(IMF) (or light cone) distributions in nucleus with parton distributions in the nucleon 
[4, 5]. On the other hand, the structure function g2(x) has no simple interpretation 
in probabilistic language [7, 8, 9] and it is natural to expect that the procedure of 
extracting g^ix) from nuclear data should be more complicated. 
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The nuclear effects for spin dependent structure functions were investigated in a 
number of papers (see, e.g. [5, 6, 12, 13] and references therein) in various model ap- 
proaches and it was shown that relativistic effects could violate the impulse approxima- 
tion. We investigate this question in the framework of Frankfurt- Strikman relativistic 
deuteron model [4]. We consider a deuterium target, but, qualitatively, our results are 
valid for any few-nucleon nuclei. 

We show that even in the IMF the qq-paiis creation or annihilation by electromag- 
netic field could violate impulse approximation for g®(x). Following the arguments 
of our previous paper [9], where in the special IMF a proper interpretation of g2(x) 
in terms of quark-parton wave functions has been proposed, we find that in such sys- 
tem the deuteron structure function g^ix) can be expressed through nucleon structure 
functions by some generalized convolution equation. We find that g^ix) is determined 
not only by g% (x) but also by derivative of xg^(x). The second term arises due to 
the nucleon internal motion and vanishes for noninteracting nucleons. The numeri- 
cal estimates show that corresponding contribution is small but increases strongly at 
X — > 1 and may achieve about 10% at X ~ 0.75. For the transverse structure function 
g^{x) the role of such term strongly depends on the twist-3 contribution and could be 
considerably larger. 



2 Kinematics. Impulse Approximation 

The electron deep-inelastic scattering on spin-one target can be described by the hadronic 
tensor: 

^-W^(\', A; P, q) = E(2vr) 4 (P, \'\J,\X)(X\J„\P, X)6(P + q - P x ), (1) 
4tt Y 

where P and q are deuteron and virtual photon 4-momenta, A and A' are components 
of initial and final deuteron spins along z-axis, J ^ is the electromagnetic current. For 
spin-one target the electromagnetic tensor depends on eight invariant structure func- 
tions [10]. We are interisted only in g\(x) and g2(x) } which are determined through 
antisymmetric part of hadronic tensor in the same form as for the spin 1/2 target: 

^W^(X', A; P, q) = iMe^ Xa q X {M 2 s°G 1 + (P • qs° - s ■ qP a )G 2 ), (2) 

47T 

where the 4- vector characterizes the target polarization state and is expressed in the 
following form 

*" = -^ a ^i(P)^(P) (3) 
through initial and final deuteron polarization vectors <p* a (P) and Lp^(P) 

P^,(P) = 0, ^(PK(P) = -1. (4) 

In terms of old fashioned perturbation theory the impulse approximation for the 
hadronic tensor (1) corresponds to the diagram of fig. la, which is supposed to dominate 
in IMF, the diagrams with g^-pairs creation or annihilation by electromagnetic current 
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(fig. lb) being suppressed at P — > oo. But this argument is not valid for g2(x). Let 
us demonstrate this statement in the conventional frame, where photon has pure z- 
component [11]: 

P, = (£,P, 0,0), q, = (0,-2Px,0,0). (5) 

In this system spin dependent SF can be expressed through antisymmetric parts of W-j 
and W°q (i = 1,2) as follows: 



— W a 

2tt tJ 



2ie ij g 1 (x] 



so 
2P 



(6) 



where the functions g\ = M 2 vG\ and g 2 = Mv 2 G 2 scale in the Bjorken limit, M is the 
deuteron mass, v = P ■ q/M } x = Q 2 /2P ■ q; we also use more conventional variable 
X = 2x related to one nucleon in the deuteron. 



P, 7 





fig. 1: Deep inelastic scattering on deuteron: a) diagram corresponding to impulse approxi- 
mation, b) diagram which could violate impulse approximation. 

For the quark momenta (which are defined on diagrams) we introduce standard 
parameterizations: 

P = riP + k 1 _ k L P = 



Pi 



XlP+PlL, p' 1 = -XlP+PlJ_, p"=-X 1 P-p 1 _ L , 



PilP = 0. 



(7) 



The vertices of photon interactions with quarks (on fig. la) and with qq-p&iis (on fig. lb) 
at P — > oo behave as follows (z,j = l,2): 

u(Pihi u (Pi) = 2Px 1 a l a 3 , u(p / 1 )7ou(pi) = 2(m + le^a.pu), 



u (Pi)n v (p") = 2(ma t + ie i jp 1 j)a 2 , uip'^ovip") = 2Px 1 a 3 a 2 . (8) 

The amplitude of antinucleon interaction with deuteron on fig. lb is ~ P. Hence, 
the large energy denominators corresponding to dashed lines on diagrams of fig. lb 
may be compensated for Wio (but not for Wij ). These diagrams may contribute to 
(gi(x)-\-g2(x)) and, therefore, may violate the impulse approximation. 
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Let us consider now IMF where photon have pure transverse component at P —> oo: 

2 2 

P, = (P,P, 0,0), 1» = (^^U)- (9) 

In this frame the structure functions g\(x) and g2(x) are expressed through antisym- 
metric components of the hadronic tensor W t a (i = 1,2) as follows: 

2Mx 

H± W? (X, A) = 2ie ij q j {g 1 (x)X8 yx + —^-g 2 (x)(s ± u) yx } (10) 



47rPx q]_ 

The structure function </2(^) can be extracted from the spin-flip amplitude, (A' = ±1, 
A = 0), So = 0, sj_q±_ = l/\/2(q x — i\q y ). In the coordinate system (9) the quark- 
photon vertices behave as follows: 

u(p' 1 )7 u(pi) = 2Pxi, u(p' 1 )7 i u(pi) = 2(pu + q, + ie lk q k a 3 ), 

u(Pihov(p") = <?<fcr 3 cr 2 , u(Pi)^ v {Pi) = -1Px 1 a i a 2 - (H) 

The contribution of fig. lb diagrams does not vanish at P — > oo. But the corresponding 
values of spin-flip amplitude for such diagrams are independent on q direction and 
therefore do not contribute to g2^x). It is easy to see if one compare terms proportional 
to q x q y at both sides of (10). For the spin-nonflip part of the hadronic tensor, which 
determines g\(x) these diagrams also do not contribute in Bjorken limit. 

Having been convinced that the diagrams with qq pair creation do not contribute 
to gP(x) and g®(x), we do not need to invoke their parton representation and derive 
the impulse approximation immediately for the structure functions. According to fig. la 
diagram, we can derive in IMF (9) the deuteron hadronic tensor through nucleon one : 

r rlV 

W^(P,A',A) = £ -^V^K^P^sl), (12) 

N,s 1 ,s' 1 ,s 2 J V 

where the sum is running over N = p } n and dT is two particle phase space: 

= , v (13) 

(2tt) 3 2t7(1 — 77) V ; 

Note, that we neglect the admixture of pions, gluons, or six quark states in the 
deuteron IMF wave function. 

The energy denominator, corresponding to the dashed line on the fig. la diagram is 
included into the definition of deuteron WF: 



- ^k'li-k) (i4 » 



where r„ „ is the deuteron-nucleon vertex function. 

The deuteron and nucleon Bjorken variables x = q\/2P ■ q } and x = q\j2Pq in IMF 



(9) are connected as follows: 

XT] 



1 + 2 -±4±i 



(15) 
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We emphasize the term 2k±q±/q 2 _ in (15) is kept in the scaling limit. Later we shall 
see that such term defines nontrivial contribution to g®(x). 

Equation (13) determines the impulse approximation. Similar equation has been ac- 
tually obtained in [5] but in terms of probabilities and therefore cannot be immediately 
used for extraction of g®(x). 

The following circumstance should be noted. The equation (13) is not an exact 
equation even in the special IMF (9). It could be violated for spin-flip amplitudes, and 
should be understood only under definite prescriptions. For the antisymmetric part of 
Woi it is valid, for instance, only in terms proportional to q x q y . 

The nucleon in our approach is off energy shell and its structure functions may 
depend on its off-shellness, i.e. on k\jm\ and on 1 — 2rj. The off-shell effects were in- 
vestigated in [12, 13]. We ignore such dependence and take into account only relativistic 
corrections which arise due to proper definition of deuteron IMF WF. 



3 The deuteron IMF wave function 

For numerical calculations the deuteron IMF wave function which is defined according to 
Eq.(14) should be connected with phenomenological rest frame WF. We do not discuss 
different approaches (see [13, 14] and references therein) and follow Frankfurt- Strikman 
prescription [4]. 

It is convinient to write the spin-orbital part of deuteron WF via Melosh matrices, 
which determine their transformations from rest frame into IMF. [15, 16]. 

K, S2 (v,kl) = J2 $( M o) U Sl , s[ ( V ,k ± )U S2 , s ,(l- V} -k ± ) ^(v.kl) (16) 



s' s' 



where $(Mq) is radial part of deuteron WF which is supposed to depend only on one 
argument- invariant mass of two-nucleon system [4] 

Ml = (17) 

where m stands for nucleon mass. The Melosh matrix has a following form: 

U(ri, & l) = , 0/ 1 (18) 

y/(m + M f]) 2 + k\ 

The spin-orbital parts of deuteron rest frame WF for S- and D-waves we write via 
Pauli spinors w s a : 

I u\ = e}w*^a % a 2 w s 2 
/ to A e}w*^a 3 a 2 w s 2 U t k 3 

Vs lS2 = 2 ( 9 ) 

where e 8 A are deuteron rest frame polarization vectors, k = (k 3 , k±), k 3 = (rj — |) M is 
the nucleon 3-momenta in the two-nucleon rest frame. 
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The normalization of deuteron WF (14) can be fixed from electric or baryon charge 
conservation: 

' * A dT = 1. (20) 

The deuteron IMF wave functions defined according to (14), (16) and (19) coincide 
with Frankfurt- Strikman WF ([4, 5]), but is written in a more compact form. To connect 
our WF (14) with existing parametrizations, we introduce conventionally normalized 
wave functions: 



U(k) = - 



L2 

7T ®u(Ml 



\ 2u 



W(k) = - 

7T 



\ 2u 



(U 2 + W 2 )dk = 1. 



(211 



4 Generalized convolution equation for g\(x) and 



To derive the deuteron structure function consider Eq.12 for antisymmetric part of W t a 
and compare terms proportional to q x q y . Making use of Eq.(10) and analogous equation 
for the nucleon hadronic tensor we find in terms of deuteron IMF wave functions (16): 



J j] 



(22) 



q —, pL 2Mg?(x) = I *VV?)* — <${*) + ** A -3*° 4 xg?(x)}, (23) 
y2 J J] J] x z 

where gf (X) = gfiX) + gf(X). The second term in (23) at first sight seems do not 
contribute to g®(x), because the value of ty* x <j 3 ty° is proportional to linear power of the 
nucleon transverse momenta ~ (k x — i\k y ) and vanishes after taking integration on k±. 
But if we take the exact equation (15) connecting deuteron and nucleon Bjorken vari- 
ables and expand xg^(x) at q\ — > oo keeping terms ~ k^q^/q 2 ^, we find nonvanishing 
contribution of this term determined by first derivatives of xg^(x). 

Note, that relativistic corrections are expected to be small (kj_/m 2 ) ~ 1/200) and 
we neglect terms proportional to k\jm 2 as compared to unity but keep them where 
they are enhanced by a factor 1/(1 — X) at X — > 1. We also keep terms proportional 
to k 3 /m. The neglected terms are essential in the range kinematically forbidden for 
free nucleon (X > 1). If necessary, they can be easily derived from Eqs.(22),(23). We 
present the simplified equations: 

9 fV) = £ ^ (*« - ^f) K*. + ^§> - S?(i) (24) 
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[-- 5-^+^* u ]^[%r(5)]} (25) 

where $„(M 2 ) and $^(M 2 ) are S- and D- wave radial wave functions. Note, that 
contribution to g^ix) arises due to D-wave part of deuteron wave function (the hrst 
term in the square brackets in (25)) and in S'-wave part due to Melosh transformation 
(second terms). Second term in square brackets in (24) is also originated from Melosh 
matrices. Eq.24 coincides with results of Ref.([5]. 

It is easy to see that for g^ix) Eq.(25) satisfies the Burkhard-Cottingham sum rule 

[17] 

[ 1 g 2 (x)dx = 0. (26) 
Jo 

Eqs.(24) and (25) can be considered as a generalization of naive convolution equa- 
tion. Due to relativistic effects the deuteron spin dependent SF cannot be presented 
in pure probabilistic form in terms of deuteron rest frame WF. Together with (x) 
contribution to g®(x), due to Melosh transformation a nondiagonal terms involving the 
interference of S- and D-wave part of the deuteron WF also contribute to deuteron SF. 
Note also an extra factor of 1/rj in (25). 

It is well known that g2(x) vanishes for noninteracting constituents. In that sense 
Eq.(25) has a transparent interpretation. The hrst term corresponds to quark gluon 
interactions forming a nucleon and vanishes for noninteracting partons, second term 
corresponds to nucleon interaction forming a deuteron and vanishes for noninteracting 
nucleons. 

Our results, in general, resemble results of Refs.[12, 13], where the binding effects 
were investigated in a covariant framework. But it is difficult to establish a detailed con- 
nection between our aproach and that of [12, 13], because underlying model assumption 
are different. 




fig. 2: Relativistic corrections to g®{X), A^X), A 2 (X), A 3 (X) and to g D L (X), A(X) (in 
%). 

In Ref. [6] the nuclear effects for g^{x) were investigated in the framework of operator 
product expansion on light cone. The authors have considered a model of nucleon 
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binding by a pion field. It is shown that impulse approximation for g^ix) is violated 
by nondiagonal contributions of states with different number of pions in the initial and 
final deuteron wave function. Keeping in mind that the space-time picture of bound 
state depends on the coordinate system, a following connection of our result with result 
of Ref. [6] can be established. In the IMF (9) a part of NN7T interaction which brought 
to nondiagonal contributions in [6] are hidden, in our approach, in the deuteron wave 
function and brought to second term in (25). The remaining part of interaction should 
be taken into account explicitly in terms of pion admixture in deuteron wave function. 
It seems that we take into account the essential part of NN interaction which can be 
approximately treated as potential, and Egs.(24) and (25) give good generalization of 
convolution equations. 

We present numerical estimates only for illustration. Relativistic effects turn out to 
be small but increases considerably at X — > 1. According to model calculations [18] 
the twist-3 contribution to g^ (x) at x > 0.5 is, possibly, small. In numerical estimates 
we take for simplicity g 1 ^ (x) ~ (1 — x) 3 and define g^ {x) through Wandzura-Wilchek 
relation [19]. We perform calculations with deuteron WF of the Paris NN potential 
[20]. The relative contributions of different terms in Eqs.(24) and (25) are presented on 
fig. 2; Ai(X) is smearing correction to g® (arises in first term of (24)), A 2 (X) is relative 
contribution of g^ (x) to g^(X), A 3 (X) is relative contribution of the second term in 
(24), which arise due to Melosh transformation. At X ~ 1 this quantities behave as 
follows: 

3P X 2 3P X 

MX) ~A 2 (X) -f^j- A 3 (X)~pj-^-. (27) 

lra z (1 — Xy \ra z 1 — X 

Thus at X ~ 1 the g^{x) contribution to g^{X) turns out of the same order 
than typical smearing correction and at lower values of X even exceeds it. It becomes 
noticeable at X ~ 0.5 (~ 2%) and achieves 10% at x ~ 0.75. 

The structure function g2(x) can be determined by measuring transverse asymmetry, 
which defines transverse structure function g±(x) = g\(x) + g 2 (:c), and the role of this 
correction is more important for g±(x). On fig. 2 we also present the relative contribution 
of second term in (25) to g^iX), A(X), arising from (xg^(x)y. If twist-3 contribution! 
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really small, g^(x), will be small as compared with g^ (x) and A(X) could considerably 
exceed A 2 (X). In our example we neglect the twist-3 part of g% (x) and A(X) turned 
out too large. Nevertheless it is clear that this correction could be very important in 
extracting g^ix) from deuteron data. 
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